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Abstract: The phase structure of QCD-like gauge theories with fermions in various 
representations is an interesting but generally analytically intractable problem. One 
way to ensure weak coupling is to define the theory in a small finite volume, in 
this case S 3 x S 1 . Genuine phase transitions can then occur in the large iV theory. 
Here, we use this technique to investigate SU (N) gauge theory with a number Nf of 
massive adjoint- valued Majorana fermions having non-thermal boundary conditions 
around S 1 . For Nf = 1 we find a line of transitions that separate the weak-coupling 
analogues of the confined and de-confined phases for which the density of eigenvalues 
of the Wilson line transform from the uniform distribution to a gapped distribution. 
However, the situation for Nf > 1 is much richer and a series of weak-coupling 
analogues of partially-confined phases appear which leave unbroken a Z p subgroup 
of the centre symmetry. In these Z p phases the eigenvalue density has p gaps and they 
are separated from the confining phase and from one-another by first order phase 
transitions. We show that for small enough mR (the mass of the fermions times 
the radius of the S 3 ) only the confined phase exists. The large N phase diagram 
is consistent with the finite N result and with other approaches based on M 3 x S 1 
calculations and lattice simulations. 
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1. Introduction 



The study of gauge theories in finite volume, and specifically S* 3 x S 1 , has been 
an interesting and fruitful one. A weak coupling regime is ensured by taking the 
size of the compact space to be small compared with the strong coupling scale, 
min[R s i , Rs'i] ^qcd- The theory is then non-trivial even at the one-loop level 
because the projection onto gauge invariant states introduces effective interactions 
between the gluons [1-4]. The large N limit is taken in order to ensure that a 
thermodynamic limit exists and genuine phase transitions occur. One motivation is 
to study thermal properties of gauge theories, in which case the S l is interpreted 
as the "thermal circle" and fermions have anti-periodic boundary conditions. A 
one-loop calculation is then sufficient to uncover the weak-coupling manifestation of 
the confinement/de-confinement transition in that on one side of the transition the 
expectation value of the Polyakov loop vanishes — the confining phase — while on the 
other side — the de-confined phase — it gains a VEV. The transition is Hagedorn-like 
in that the density of states grows exponentially in the low temperature phase. In 
order to ascertain the order of the transition higher loop effects are crucial [2]. In 
pure gauge theory it is known to be a first order transition that occurs at a lower 
temperature than the Hagedorn transition in the non-interacting theory. The same 
kind of transition occurs when adjoint matter fields are added, and in particular for 
the M = 4 gauge theory. One of the deep insights to emerge from the AdS/CFT 
correspondence is that this transition can also be seen in the strong coupling gravity 
dual as a Hawking-Page transition from thermal AdS space to an AdS black hole [5]. 
In the present paper, since we are interested in theories with periodic boundary 
conditions for fermions around the S 1 the phase transitions are quantum rather than 
thermal and the connection with the AdS/CFT correspondence is not so obvious 
even though one could imagine obtaining our theories from an M = 4 theory with 
SUSY breaking mass deformations. 

The phase diagram of SU(N) gauge theories in finite volume can be studied 
in other ways. At strong coupling lattice simulations are the dominant technique 
for obtaining the phase diagram. The phase diagrams using the strong and weak 
coupling techniques have so far not been easy to compare, since the relationship be- 
tween lattice bare parameters and continuum renormalized parameters is not clear. 
However, in some cases a qualitative comparison is possible. In this paper we ex- 
plore this possibility in adjoint QCD, that is, SU(N) gauge theory with fermions 
in the adjoint representation. The phase diagram of this theory is quite rich when 
considering Nf > 1 Majorana flavours with fermions of finite mass to which periodic 
boundary conditions have been applied in the temporal direction [6-12]. The intu- 
ition behind this is that the gauge field and fermion terms in the effective potential 
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have opposite signs and compete to dominate the Polykov loop, with the fermions 
having a disordering effect. For a sufficient number of light fermions, the disordering 
effect dominantes and a confining phase results with vanishing Polakov loop, how- 
ever as the masses are increased the disordering effect becomes weaker and a phase 
transition can occur where the centre symmetry is (partially) broken. 

The issues that we investigate in this paper are relevant to some active areas of 
research. The application of periodic boundary conditions to adjoint representation 
fermions causes the confined phase to be accessible at weak coupling [6,8,13] l . This 
is interesting for several reasons. For one, it is possible that observables may not 
differ significantly in the pure Yang-Mills theory confined phase and the perturbative 
confined phase of adjoint QCD. The latter feature might naively be inferred from 
the results of lattice calculations of the string tensions in the confining phase of 
the pure Yang-Mills theory which show little temperature dependence. However, it 
is also known that the string tensions have a significant temperature dependence 
above T c in the de-confined phase [14]. This suggests an important question: Is the 
temperature dependence a result of the change in confinement scale when moving 
above T c such that it would also occur in a high temperature confined phase, or is it 
something intrinsic to the de-confined phase? In the case where the high temperature 
confined phase is induced by periodic boundary conditions the question turns into 
one of dependence of observables on the length L of S* 1 in the weak-coupling confined 
phase. In this case the question of temperature dependence becomes one of volume 
dependence. 

The idea of volume dependence — or rather independence — is particularly inter- 
esting in the context of large N confining theories. In 1982 Eguchi and Kawai pro- 
posed volume independence in large N Yang-Mills theory in [15] where they employed 
large N factorization to show that pure Yang-Mills theory formulated on a lattice 
at some arbitrary volume, can be mapped onto the theory formulated on a single 
site. Around that time it was also shown that volume independence can only hold 
if certain symmetries are not broken, in particular it can only hold in the confining 
phase of large N gauge theories [16,18-20]. Several ideas were proposed to maintain 
the Z^r symmetry which is required on a space with d independently compactified di- 
mensions (i.e. where certain dimensions R are compactified to S l ). Two of the most 
promising proposals were the quenched [20], and twisted [22] Eguchi-Kawai models. 



x It should be noted that even though the confined phase is perturbatively accessible in adjoint 
QCD when min[i?si , R s s] <c A -1 , it was found in [33] that semiclassical analyses on M 3 x S 1 , 
specifically including the contribution of magnetically charged objects in the confined phase, are 
only valid when RgiNA <C 1. For certain observables which exhibit volume independence in the 
confined phase this should also be true on S 3 x S 1 . We thank Mithat Unsal for pointing out this 
important result. 
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In the quenched EK model the eigenvalues of the Polyakov loop quenched, that is 
their eigenvalues are set by hand, such that 7i% symmetry required is explicit. How- 
ever, results from lattice simulations of SU(N) gauge theory in [21] show that large 
N reduction using the quenched EK model breaks down. In the twisted EK model 
the boundary conditions are twisted by multiplying each plaquette in the lattice ac- 
tion by an element of the center such that the action becomes invariant under the 
required symmetry. The twisted EK model was shown to break down in [23,24]. 
More recently, the idea of large N volume independence was picked up again and 
generalized in [13] where the authors proposed that since QCD(Adj) has a confined 
phase which is perturbatively accessible, a generalized Eguchi-Kawai large N re- 
duction could relate the weak-coupling, small volume confined phase, to the strong 
coupling, large volume confined phase. Since this proposal there have been several 
tests. In [9] the authors performed a weak-coupling calculation of the effective po- 
tential for a three-dimensional adjoint matter theory which is related by generalized 
orbifold projection to QCD(Adj) in four dimensions with one dimension compact. 
In [11] it was shown that this calculation is renormalization scheme dependent and 
that for Eguchi-Kawai reduction to hold it requires double trace counter-terms with 
coefficients defined so that the Zjy symmetry is preserved in the limit of zero adjoint 
fermion mass. In a lattice simulation [10] for N = 3 the authors calculated the phase 
diagram of QCD(Adj) on a 16 3 x L c lattice and showed that there is a confined phase 
at both strong and weak coupling but that they are separated by phase transitions 
into other phases depending on the value of the adjoint fermion mass. From their 
results the presence, or lack thereof, of phase transitions in the chiral limit is unclear. 
But, their result places boundaries on the validity of large N reduction either way 
since it can only hold when the theories to be mapped are in confining phases. Most 
recently in [12] the Eguchi-Kawai reduced (single-site) model with dynamical adjoint 
fermions was studied on the lattice and it was shown that the Z^ symmetry is un- 
broken for light enough fermion mass. Some applications of Eguchi-Kawai reduction 
can be found in [25, 26] where the authors use large N reduction of SYM (J\f = 1 
in [25] and Af = 4 in [26]) on S 3 x 1 to reduce the theory to a single dimension for 
the purpose of studying supersymmetric matrix quantum mechanics. 

Our calculations in this paper place boundaries on validity of large- N reduction 
at weak coupling and finite fermion mass by mapping the regions of Z(N) symme- 
try breaking in the phase diagram of large iV adjoint QCD. The phase diagram is 
calculated as a function of the length L of S 1 , the radius R of S* 3 , and the adjoint 
fermion mass m. In particular, for small L/R the mass of the fermions, times the 
radius R, must be below a critical value to keep the theory in a Z(N) symmetric 
phase. It is clearly an important question to understand the phases of QCD with 
adjoint fermions as a function of the volume and mass and it is to this question that 
we now turn. 
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In Section 2 we compute the effective action for the theory on S 3 x S 1 as a 
function of the Polyakov loop to one-loop order paying particular attention to the 
inclusion of a mass for the fermions. Section 3 investigates the phase diagram of the 
large N theory as a function of the radii of S 3 x S* 1 and the mass of the fermions 
and for different numbers of adjoint fermion flavours. Here we show the existence 
of a rich phase structure for Nf > 1. In the final section, we consider the same 
theories with N finite where the phase transitions are no longer non-analytic but 
allow comparisons with with earlier work mentioned above for the theory on 1R 3 x S 1 
and lattice simulations. 



2. The Effective Action on S 3 x S 



In this section, we review the way that the effective action is calculated on S* 3 x S 1 to 
the level of the one-loop approximation. The only new ingredient over earlier work is 
the inclusion of a mass for the fermions. Our approach follows closely the philosophy 
set out in the beautiful paper [1], however, we shall use a more conventional form of 
gauge fixing, described in [27-29], that leads to the same result. 

We shall start with SU (N) gauge theory with a number of Majorana fermions 
tpf transforming in representations Rf of the gauge group. 2 The action is 

i r Nf 

S = — rf 4 x v ^{-iTrF^ + ^(#/^/-m / ^ / )} (2.1) 

and the covariant derivatives are appropriate to the representation Rf. 

The problem before us is to compute a Wilsonian effective action for the gauge 
theory on S 3 x S* 1 to the one loop order. We denote the length of S* 1 by L and the 
radius of S 3 by R. The only zero modes in the set-up belong to A , the gauge field 
component around S 1 : 

a = ^ [ A . (2.2) 

VolS 3 x S 1 J s t xS i 

We can use a global gauge transformation to diagonalize a: 

a = L -1 diag(0i) . (2.3) 

The 9i are angular variables since there are large gauge transformations (but periodic 
around S* 1 ) that take 6i — > Qi + 2n. Physically, the gauge invariant quantities are 

2 In general in order to have a mass term the representations Rf must be real or include complex 
conjugate pairs. In the case here we are considering the adjoint representation which is real. 
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built out of the Wilson loop P = e tLa (the Polyakov loop in the case of thermal 
boundary conditions) evaluated in the fundamental representation: 3 

N 

Tr P n = eind] , ( 2 - 4 ) 
i=i 

On top of this there are additional large gauge transformations that are only periodic 
on S 1 up to a subgroup of the centre Z N depending on the matter content of the 
theory. For adjoint matter, this subgroup is the whole of 1*n and so non-periodic 
large gauge transformations take 9{ — > Q{ + 2n/N and so transform TrP by an N- 
th root of unity. Hence, strictly speaking, the gauge invariant observables are, for 
example, |TrP|. 

The radiative corrections at the one loop level are obtained by taking the constant 
mode ( |2.3[ ) as a background VEV and integrating out all the massive modes of the 
fields. To this end, we shift A — > A + a and then the one-loop contribution involves 
the logarithm of the resulting functional determinants which depend on a in a non- 
trivial way. As usual we have to fix the gauge in some way and we prefer to use a 
different and more conventional approach than that of [1]. To this end we impose 
Feynman gauge by adding the standard gauge fixing term 

Sgf = 4~\ [ d'x^Trh^ + D A )\ (2.5) 
9ym * J 

and appropriate ghosts. To the one loop level, we expand the action to quadratic 
order in fluctuations. The gauge field part of the action, including the ghosts, is 

Sgauge = ^-[ d 4 x^Tr \\A,{-Dl - A^)A 

9ym J 1 ^ (2.6) 

+ ~ A(V) ) A ' + +c(-D 2 - A< S V 

Here, A^ and A^ are the Laplacians on S* 3 for scalar and gauge fields, respectively. 
The scalar Laplacian is simply A^ s ^ = g -1 / 2 ^^ 1 / 2 ^) whilst the vector Laplacian is 

A (v) A i = Vi V j ^ - KjA? , (2.7) 

where Rij is the Ricci tensor of S 3 . In the above, D = do + iot and so includes the 
coupling to the VEV 

Each fluctuating field is expanded in terms of appropriate harmonics on S 3 x S 1 
and a typical contribution to the effective action is of the form 

±Tr R \og(-D%- A) , (2.8) 



3 Unlcss otherwise specified, traces are taken in the fundamental representation. 
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field 


angular mom. 


energy 


degeneracy 




e > o 




2£(£ + 2) 


Q 


£ > 


y/£(£ + 2)/R 


(£+l) 2 


c, c 


£ > 


[/£{£ + 2) /R 


-2(£+l) 2 


A 


£ > 


y / £(£ + 2)/R 


(£+l) 2 




£ > 


^J(£ + \) 2 + m 2 R 2 /R 


-2£{£ + l) 



Table 1: The fields, their angular momenta, energy and degeneracy (with ± sign for c- 
number and Grassmann fluctuations) in the effective action. The fermion result is for a 
massive Majorana fermion on S 3 x S 1 . 



the ±1 being for c- number and Grassmann fluctuations, respectively and A is the 
Laplacian on S 3 appropriate to the tensorial nature of the field on S* 3 ; either A^ s \ 
/\(v) or The background VEV a acts as a generator of the Lie algebra of 

SU(N) in the representation R of the gauge group appropriate to the field and the 
trace includes a trace over that representation of the gauge group. The eigenvalues 
of do are simply 2mn/L, n G Z, while the eigenvectors of the Laplacian on S* 3 are 
labelled by the angular momentum £: 

Aij e = -e 2 e ij e , (2.9) 

and we denote their degeneracy as de- The ee and di depend on the field type. We 
review the spectra of the appropriate Laplacians on an arbitrary sphere in Appendix 
A. For us, the relevant fields are scalars (more precisely minimally coupled scalars), 
vectors and spinors and below we list the relevant data: 

(i) Scalars. For minimally coupled scalars e e = R~ l y '£{£ + 2) and the degener- 
acy d t = {£+l) 2 with £ > 0. 

(ii) Spinors. For the irreducible 2-component real spinors, 4 we have Ei = R~ 1 (£+ 
1/2) and d e = £(£ + 1) with £ > 0. 

(iii) Vectors. Here the situation is more complicated. A vector field V* can be 
decomposed into the image and the kernel of the covariant derivative: V* = ViX + Bi, 
with V l Bi = 0. The eigenvectors for the "transverse part", B i: have Et = R~ l (£ + 1) 
and d t = 2£(£ + 2) with £ > 0. On the other hand, the "longitudinal part" V^x has 
e t = R~ 1 y / £(£ + 2) with degeneracy d e = {£ + l) 2 but with £ > only. 

It is a standard calculation using the identity f|^° =1 (l + x 2 /n 2 ) = sinh(7ra;)/(7ra;) 
4 A Majorana spinor on S 3 x S 1 corresponds to 2 such spinors on S 3 . 
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to show that ( j2.8|) is equal, up to an infinite additive constant, to 

oo oo 1 

^^|L^-2^-e-^Tr R (P")} . (2.10) 

1=0 71=1 

The first term here involves the Casimir energy and since it is independent of a will 
play no role in our story and we will subsequently drop it. 

Notice that the t > contributions from Aq, Ci and the ghosts all cancel leaving 
only a net contribution from the I = modes of the form 5 



oo _ 

E-^( pn ) • ( 2 - 12 ) 



n 

71=1 



The remaining gauge modes are the vector modes E>i and the fermions. Using the 
sum (|2.10|) and including the Jacobian term in ( |2.12|) , the full effective action is then 



oo Nj 

S ^ = J2-{( 1 ~ Zv(nL/R))Tr^(P n ) + J2z f (nL/R,m f R)Tr Rf (P n )j . (2.13) 

71=1 H f=l 

In the above, we have defined 

z v (L/R) = 2£/(/ + 2) e -*»V = fe (1 _ e J; B)3 , (2.14) 
and for the spinors 

oo 

z f {L/R, mR) = 2j2 W + l) e -V^+V2) a +m^/fl (2 . 15 ) 
i=i 

We have assumed that the fermions have periodic boundary conditions around S . If 
one wanted to describe the case of finite temperature the fermions have anti-periodic 
boundary conditions and Zf(nL/R,m,fR) must be multiplied by an additional factor 
of (-l) n . 

Notice for the vector modes, we are able to perform the sum over the angular 
momentum, but for the fermion modes this is not possible due to the non-zero mass. 



5 This part is precisely the exponentiation of the Jacobian that converts the integrals over the 9i 
into an integral over the unitary matrix P = diag(e i ): 

„ N oo_. „ N ft ft f 

/ [] Mi exp { £ ^Tr adj (P«) } oc / J[ d9 t J[ sin 2 (^^) = / dP . (2.11) 

i=l n=l i=l i<j 

However, we will leave the Jacobian in the exponent since it must be considered as part of the 
effective action for the eigenvalues. 
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However, a useful expression for the fermionic contribution can be obtained by apply- 
ing a version of the Abel-Plana formula which is proved in Appendix B appropriated 
to a function with branch points on the imaginary axis: 

e/(*+d= r dx M - i T dx f{ix + y * ~ g) ■ ( 2 - ie ) 

In the above e is positive, real, and infinitesimal. Applying this formula to the 
function 

/(f) = 2t{t + 1 /2) e - L ^^ +1 W 2+m2R2 l R , (2.17) 
gives the integral representation 
Zf(L/R, mR) 

poo poo ™2 _|_ 1 

= 2 dx(x 2 - k) e -Wx2+m*&/R + 4 / rfa; _j: 4 si^Ly/x 2 -m?R?lR) 

JO JmR 6 nx + 1 

m v> r°° x 2 + •! , 

-K 2 (Lm) - —K x (Lm) + A J dx — + ^ sin(L v / x 2 - m 2 R 2 /R) . 



L 



mR 

(2.18) 



Note that in our case the function /(ix) is real for x < mi? and becomes imaginary 
for x > mR and so the lower limit of the second integral has been taken to be mR. 



There are two consistency checks we can make on the integral expression ( 2.18 ). 
Firstly, in the massless limit, m — > 0, we have 

A p -3L/2R °° 

z f (L/R, 0) = (1 _ e _ L/R) , = £ 2£ ( £ + l)e~ L{i+1/2)/R • (2.19) 

Secondly in the limit R — > oo with fixed m and L, we have 

2m 2 R? 

Zf(L/R, mR) — > — - — K 2 {Lm) , (2.20) 

Lj 

which is the expression that one obtains by working directly on M 3 x S 1 [8]. 



3. Phase Structure of QCD(Adj) 



In this section, we apply our result ( [2.13] ) to the particular case of a theory with N 



f 



adjoint fermions with equal masses. In this case, we have 



N 



Tr adj (P") =^co8(n(^-^)) (3.1) 
ij=l 
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and so 

oo 1 N 

S(6i) =J2~( 1 ~ z v( nL /R) + N f z f(nL/R, mjR)) ^ cos(n(0; - d 3 )) . (3.2) 

77 V / 
n=l 13=1 

The phase structure is determined by minimizing this with respect to the Since 
N is large, it is more appropriate to describe the configuration in terms of a density 
p(9) normalized so that 

dOp{0) = 1 . (3.3) 







In this case, we can write the effective action as 



s[p(6)} = n 2 f dd [ dd'p{6)p{d') tt nL l R ' mK ) cos ( n (# _ e >)) (3.4) 

J J n =l n 

In the above, we have defined the function 

f(L/R, mR) = 1 - z B (L/R) + N f z F (L/R, mR) . (3.5) 



It is useful to Fourier analyze the density: 

_ 00 



2tt 

n=— 00 



with po — 1 an d Pn = P- n for reality. The action then becomes 



q\ law N 2 ^ f(nL/R,mR) l 2 

S[p{0)\ = — 2^ \Pn\ (3-7) 

n=l 



The phase structure hinges on properties of the function f(L/R, mR) and specifi- 
cally on its sign. For large L/R and any 777, both zs,zp — > and so f(L/R, mR) — ^ 1. 
The behaviour in the limit of small L/R depends on how m is scaled. If we keep mL 
fixed then 

f(L/R, mR) — (N f K 2 (mL) - ^|_) . (3.8) 

In this limit, / is positive (negative) for mL < a (mL > a), where a is the solution 
of 

N f a 2 K 2 (a) = 2 , (3.9) 

which is only possible if Nf > 1. For Nf — 1, / is negative for all mL (in the limit 
of small L/R). If we consider small L/R but keep mR fixed then 

4R 3 m 2 R 3 Nf 
f(L/R, mR) — -g-fty - 1) - — ^ + ■ ■ ■ . (3.10) 
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Note that this is positive for Nf > 1. 

It will be useful to chart the phase diagram initially in the (L/R, mL) plane and 
in Figure [I] we show the corresponding regions for which / is positive and negative, 
for N f = 1,2,3,4. 

8 
7 
6 
5 

mL 4 
3 
2 
1 


0.3 0.6 0.9 1.2 1.5 

L/R 

Figure 1: The lines where / = indicating the regions where / > and / < in the 
(L/R, mL) plane for Nf = 1, ... ,4. 

Having charted the region where / is positive and negative we can now build up 
a picture of the phase structure. It will be useful to define f n = f(nL/R,mR) so 
that 

N 2 °° f 

s =t£t> |2 ' (3 ' n) 

71=1 

In the region where all the f n > 0, n — 1, 2, . . ., it is clear that the action will be 
minimized when all the Fourier modes vanish p n = 0, except po = 1, i-e. p(0) = 1/2-7T. 
This phase is the weak-coupling manifestation of the confining phase where the centre 
symmetry is unbroken and TrP n = 0, n > 0. In view of this we will call it the 
"confining phase" . 

The confining phase covers the region of the phase diagram where all the f p > 0. 
This is separated from the rest of the phase diagram by the union of arcs where a 
given f p vanishes. The arc where f p = extends between 2 multi-critical points 
where f p = f p+ i = and / p _i = f p = 0, respectively. As one crosses the contour on 
which / p = 0a phase transition occurs where the density p{6) develops p gaps, that 
is p intervals around the circle on which p(6) vanishes. We call the resulting phase 
the Z p phase since at strong coupling it would be identified with a partially confined 
phase where the SU(N) gauge group is broken to a subgroup SU(p) that confines 
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and a Z p subgroup of the centre symmetry remains unbroken. The signature of the 
Z p phase is that the order parameters behave as 

TrP n = 0, n/p&Z, 

/yT (3.12) 
TrP n ^0, n/peZ. 

The detailed argument of why such a transition occurs follows as a generalization 
of the transition from the uniform to one gap phase described in [1] and we include 
it for completeness. The important point is that the configuration space {p n } has 
a non-trivial boundary which encloses the allowed region surrounding the origin 
because the density p{6) cannot be negative. The boundary region consequently 
consists of distributions for which p{9) vanishes at a subset of points (including 
finite intervals) around the circle. Furthermore, it is clear that allowed region in the 
configuration space is a convex region, since if pi(0) and p2{6) are allowed then so is 
t Pl (0) + (1 -t)p 2 (6), < t < 1. 

Let us consider the transition across the line where f p — 0. On the confining 
phase side of the transition f p > and p(6) = 1/2tc. At the transition point, the 
action becomes independent of the Fourier component p p and so the one complex 
parameter family of densities 

p(e) = ^(l+p p ^ e + P ;e-^) (3.13) 

for < \p p \ < | all have the same (vanishing) action. As f p becomes negative, 
then the points with S — in the configuration space is a cone whose angle opens 
as f p becomes more negative. The locus of configurations with S < correspond 
to hyperbolae lying inside this cone. It follows that the configuration with minimal 
action lies on the boundary of the configuration space at a point where one of the 
hyperbola lies tangent to the boundary, i. e. on a distribution where p{9) vanishes at a 
subset of points. This is illustrated in Figure [| We conclude that new configuration 
is continuously connected to the \p p \ = | density in ( |3.13| ) and so must have precisely 
p gaps. By symmetry this phase will be invariant under a Z p subgroup of the centre 
symmetry. A schematic view of the transition appears in Figure |3|. The transition 
is first order since the effective action is discontinuous across the transition. In the 
confining phase we have S = 0, whereas just above the transition to leading order 
it is sufficient to take the density to be (|3.13| ) with \p p \ = | and plug this into the 
action to get 

N 2 

S = —f p (L/R,mR) , (3.14) 
4p 

because the density itself only changes at a higher order. Since f p (L/R, mR) has non- 
vanishing first derivatives in 5L and 5m 2 the derivatives of S change discontinuously 
across the transition implying that it is first order. 
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/ boundary 



Figure 2: The structure of configuration space showing \p p \ and one additional direction. 
The boundary is indicated by the dotted line and it is important that the allowed region 
is convex. For f p < 0, the lines of vanishing action define a cone and the lines of constant 
negative action being hyperbolae therein. The density (i) is the uniform distribution char- 
acteristic of the confining phase; (ii) is the density with \p p \ = ^ which lies at the boundary 
of configuration space; and (iii) is the density with minimal action as lying at the boundary 
of configuration space where the lines of constant S lie tangent to boundary. 




(i) (ii) (iii) 

Figure 3: The behaviour of the density across the transition at f p = 0: (i) the uniform 
density in the confining phase and (iii) the Z p phase (for p = 3). At the transition point 
(ii) the mode p p becomes massless and the density develops p zeros as shown in the middle. 



As we remarked above the critical lines f p = and f p+ \ = cross at multi-critical 
points where f p = f p+ i = 0, and at these points the confining, Z p and Z p+1 phases 
are all continuously connected via the density 

p{9) = i-(l + p p e^ + p* p e-^ + p p+1 e^ e + p^e'^ 6 ) . (3.15) 

By continuity, it must be that there are lines of first order phase transitions that 
separate the partially confined phases and which end on the critical points that lie 
somewhere between the continuation of the f p = and f p+1 = lines. 6 The actual 

6 The argument is as follows, in the vicinity of the critical point we have S = N 2 f p /4p and 
N 2 f p+ i/A(p + 1), respectively, in the Z p and Z p+1 phases. A first order transition occurs when 
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positions of these lines of transition will depend in detail on the gapped distributions. 
However, in the limit R — > oo with m and L fixed, we know from (|3.8| ) that up to 
an overall factor f(L/R,mR) oc g(mL). Hence, these lines of first order transitions 
must asymptote to lines of constant mL. 

In the Z p phase, as L/R — > we expect the gaps will grow and in this limit the 
density p{6) in the Z p phase will only have support at p equally spaced points around 
the circle 

1 

— - - 27Tj/i>) ■ (3-16) 

P 3=0 

In this case, the line of transitions will occur for mL being the solution of the equation 



f(jpL/R, mR) f(j(p + 1)L/R, mR) 

2^ Jp 2^ j( P + i) °' ' 



or more concretely since L/R — > we can use Q3.8Q to get the conditions 

L) 2 J 

(3.18) 



OO j 

= E f? - fp+ i^ (^((p + l)jmL)- 



-(j(p+l)) 2 V J zvv " ;J y ((p+l)jmL) 2 , 



The phase diagram in the (L/R,mL) plane is shown in Figure |] for the two 
distinct cases Nf = 1 and Nf = 2. In the former case, only the phase with one 
gap appears which at strong coupling is identified with the de-confined phase since 
the centre symmetry is completely broken. For Nf > 1 all the Z p phases appear for 
p = l,2,.... 

The transition between the confined phase and fc-gap phases occurs along the 
fk = curve. The line of transitions between fc-gap and k + 1-gap phases only 
slightly differs from the fj. = curves of Figure |j. From the M 3 x S 1 result in [8] 
the transition between the 1 and 2-gap phases occurs for mL ~ 2.020, where as 
the mL asymptote of the /i = curve lies at mL ~ 2.027. This difference is only 
barely visible in Figure |j. The confined phase is always favoured below a critical line 
with slope (mR) c which increases with Nf. The (L/R) c ~ 1.317 line in both figures 
indicates the value of the deconfinement temperature of the pure SU(oo) Yang-Mills 
theory determined in [1]. 

fp/p = fp+i/(P + 1) which must necessarily be in a region where f p < and / p+1 < 0, in other 
words somewhere between the continuations of the f p = and f p +\ = lines. 
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0.3 0.6 0.9 1.2 1.5 0.3 0.6 0.9 1.2 1.5 



(a) L/R ( b ) L/R 

Figure 4: The phase diagrams in (L/R,mL) coordinates for (a) Nf = 1 and (b) Nf = 2. 
We have shown the transitions between the Z p and phases along the continuation of 
the f p = line for simplicity since this seems to be a good approximation and matches the 
value calculated in the M 3 x S 1 from [8] indicated by the arrow pointing to the rnL-axis. 

In Figure |5|, we re-plot the same phase diagrams in the (L/R, mR) plane. Again 
it is clear that only the confining phase exists for small enough mR. The basic form 




(a) L I R (b) LjR 

Figure 5: The phase diagram in (L/R,mR) coordinates for (a) Nf = 1 and (b) Nf = 
2. As above, we have shown the transitions between the Z p and phases along the 

continuation of the f p = which seems to be a good approximation. 

of the phase diagram can be understood intuitively as follows. The important point 
is that periodic fermions contribute positively to f(L/R, mR) and consequently tend 
to act so as to disorder the Polyakov loop, counteracting the effect of the gauge field, 
and favour the confined phase. However, their effect goes away as m increases due to 
decoupling. Consequently, for large fermion mass m, the matter fields decouple and 
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we recover the confinement /de- confinement of the pure gauge theory. Whereas for 
small mass the fermions win in the competition with the gauge fields and disorder 
the Polyakov loop. The most striking result of our analysis is that the confining 
phase extends all the way down to small L/R as long as the fermion mass in units 
of l/R is below a critical value, m < 1.225/ R for N f = 1, m < 3.203/i? for N f = 2. 



4. Finite N 

It is useful to also consider QCD(Adj) at finite iV for several reasons: (i) we can 
determine how the large N limit is approached and develop some intuition of when 
finite N results start to approximate those in the infinite N limit, (ii) with finite N 
and finite volume better qualitative comparison with lattice results is possible, and 
(iii) it is possible to compare with finite N results on M 3 x S 1 by considering the 
limit R — > oo. 

To build intuition on how the infinite N limit is approached it is helpful to 
remove fermions for the moment and consider pure Yang-Mills theory at finite N. 
In the strong coupling limit the de-confinement phase transition has been observed 
in lattice simulations and it is believed that the transition is second order for N = 2 
and first order for N > 3 [30]. 

The weak-coupling analogue of the de-confinement transition of pure Yang-Mills 
theory is observable from perturbation theory on small volume manifolds at weak 
coupling. In [1] the authors calculated the de- confinement temperature from one- 
loop perturbation theory on S 3 x S 1 and found it to occur at T d R = 0.75932 in the 
large N limit. In [2] the same authors computed higher loop corrections to show that 
the large N de-confinement transition is first order in the weak coupling limit. 

By numerically minimizing the Yang-Mills effective potential with respect to the 
Polyakov loop eigenvalues it is possible to compute the trace of the Polyakov loop 
order parameter as a function of the temperature as shown in Figure [] (Left) for 
N = 2, 3, 4 (See Appendix |C| for a short discussion of numerical minimization of the 
effective potential for finite N). The discontinuity in the trace of the Polyakov loop 
is a clear indication of the de- confinement transition even for N = 2. Increasing N 
causes T^R to approach the large N result from [1] . Near to the transition the points 
are separated by AL/R = 0.01. With this resolution the N = 2 transition appears 
much smoother than for N = 3,4, as we might expect from lattice results. However, 
it is not possible to decipher the order of the transitions for certain without taking 
the infinite volume limit. 
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Figure 6: Yang-Mills theory for N = 2,3,4: (Left) (TR, TrP). For N = 2: 0.549 < 
T d R < 0.552, N = 3: 0.641 < T d R < 0.645, N = 4: 0.676 < T d i? < 0.680; (Right) (TR, 
Vym)'- Vym is the free energy density minus the const/TP Casimir term. For larger values 
of N this result appears increasingly compatible with a first order transition. 



The corresponding Yang-Mills effective potential (eq. (|3.2|) with Nf = or 
m = oo, which doesn't include the const/Ti? Casimir term) is shown in Figure 
|6| (Right). The slight hump (most visible in the N = 4 result) indicates the ap- 
proximate location of the de-confinement transition. It appears to become more 
well-defined with increasing N. Even though the finite iV transitions shown in fig- 
ure |B| are not genuine phase transitions since this system does not have an infinite 
number of degrees-of-freedom, true transitions can be obtained by taking iV or R to 
infinity. The sharpness of the hump in the effective potential serves as an indicator 
of how well the finite N transition approximates a true infinite N, or infinite volume 
transition. 

It is important to emphasize that the use of the saddle point approximation to 
determine the preferred configuration of the gauge field is not strictly valid in the 
limit of finite N. However, one can show that it is still a reasonable approximation 
for purpose of obtaining the phase diagram, even for iV = 2, 3, by plotting the ef- 
fective action in the complete configuration space of the Q{. The effective action has 
clear minima in the configuration space corresponding to the eigenvalues obtained 
with the saddle point approximation. The competition with other configurations is 
minimal. For L/R below (L/R) c there is a clear distinction between the confined and 
gapped phases which can be shown by performing the integrals over the gauge fields 
to obtain the partition function, Z, then plotting e _5 Tr P/Z in the full configuration 
space and considering a small radius around the values of the 6i corresponding to the 
minima of the effective action. Above (L/R) c the fluctuations of the Polyakov loop 
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0.5 1 1.5 

L/R 

Figure 7: Phase diagram of QCD(Adj) for N = 3, N f = 4: (Left) in the (L/R, mR) 
plane with 3.6 < (mR) c < 3.8; (Right) in the (L/R, mL) plane with 2.5 < (mL) c < 2.75. 

in the full configuration space of the 0j are more severe, but taking the average in a 
small radius around the configurations determined by the saddle point approxima- 
tion still results in (TrP) = 0. The fact that the integrals over the gauge field can 
be solved numerically serves as a check of the results of the saddle point approxima- 
tion. However, expectation values of the Polyakov loop, for example, will always be 
zero, so it is important to find the eigenvalues from the saddle point approximation, 
and compare with plots of e~ s TrP/ J[d6]e~ s , as a function of the (9j, such that the 
different Z(N) phases can be distinguished. Using the saddle point approximation 
serves as a means of obtaining sharper phase transitions than would be observed 
by performing the full integrals over the gauge field configurations. It additionally 
picks out eigenvalues for a single minimum of the effective potential, avoiding the 
issue of finding a suitable order parameter for distinguishing the Z(N) vacua (or the 
relevant subgroup) which would otherwise be necessary since averaging over the full 
configuration space gives (TrP) = in all phases. 

For QCD(Adj) at finite N we can perform an analysis similar to what was done 
in the case of large N. The effective potential has the exact same form except that we 
can't solve the path integral using the saddle point approximation unless N is large 
enough. Since we don't know when this is true we again use numierical minimization. 

Consider SU(3) QCD(Adj). The phase diagram for Nf = 4 Majorana flavours 
in the (L/R,mR) plane is shown in Figure [I] (Left), and in the (L/R,mL) plane in 
Figure [7] (Right). The phases are defined according to the value of TrP. The confined 
(red) phase has TrP fc = for all k. The 2-gap (green) phase can be distinguished 
from the de-confined or 1-gap (blue) phase in either of two ways: 1) The 2-gap phase 
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Figure 8: Phase diagram of QCD(Adj) for N f = 2 in the (L/R, mR) plane: (Top Left) 
N = 3: 2.6 < (mR) c ^ d < 2.8; (Top Right) iV = 4: 2.8 < {mR) c ^ d < 3.0; (Bottom Left) 
N = 5: 2.8 < [mR) c ^ d < 3.0; (Bottom Right) iV = 6: 3.0 < (mR) c ^ d < 3.2. The f n = 
curves indicate the lines of transition for the N = oo result. 



has |TrP 2 | > |TrP|, whereas the de-confined phase has |TrP| > |TrP 2 |, 2) The 2-gap 
phase has Proj Z3 TrP < whereas the de-confined phase has Proj Zg TrP > 0, where 
Proj Z3 indicates projection onto the nearest Z3 axis. Darker shading in a A;-gap phase 
indicates a greater mag nitude of |^TrP fc |. The pure Yang-Mills theory transition 
is visible at large mR for 1.55 < (L/R) c ^ d < 1.6 in approximate agreement with 
the N = 3 result in Figure || For small enough mR phase transitions are avoidable 
for all L/R. The critical value (mR) CT \ t below which symmetry-breaking phase 
transitions are absent is less for smaller Nf, but the overall shape of the phase 
diagram is otherwise qualitatively similar. 

To see how the large N phase diagram unfolds it is useful to also consider the 
phase diagrams for N = 4,5, 6. To compare with the 1R 3 x S 1 results in [8] we consider 
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Figure 9: Phase diagram of QCD(Adj) for Nf = 2 in the (L/R, mL) plane: (Top Left) 
N = 3: 1.5 < (mL) c < 1.75; (Top Right) N = 4: 1.0 < (mL) c < 1.25; (Bottom Left) 
N = 5: 0.75 < (mL) c < 1.0; (Bottom Right) N = 6: 0.5 < (mi) c < 0.75. The /„ = 
curves indicate the lines of transition for the N = oo result. 



Nf = 2. The phase diagrams of iV = 3, 4, 5, 6 QCD(Adj) are shown in the (L/R, mR) 
plane in Figure |8|. As A" is increased a new phase is formed for each odd N. For N = 3 
this phase is the 2-gap phase with |TrP 2 | > |TrP|. For A^ = 5 the new phase is the 
3-gap (yellow) phase with |TrP 3 | > |TrP 2 | , |TrP|. For even A^ the phases fan out 
into the region of small L/R. But, for N > 7 the phase diagram gets even more 
complicated. For example, from [8] we know that for A^ = 7 it is possible to have two 
different 3-gap phases, one which maximizes |TrP 3 | and the other which maximizes 
|TrP 4 |. However, if JV mod k = then there is always a k-gap phase for which the 
eigenvalues are distributed in k evenly spaced clumps containing N/k eigenvalues 
per clump such that TrP fc ^ and TrP' = for I ^ k. UN mod k ^ then there 
is still a k-gap phase with k clumps of eigenvalues and for which |TrP fc | > |TrP'| 
for I not a multiple of k. We suspect the other types of fc-gap phases to result as a 
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consequence of the decreased symmetry of the finite N theory, and conjecture that 
they should not be present in limit iV — > oo. 

As iV is increased the new phases extend down into lower values of (mR) than 
the de-confined phase, however when the next phase is formed the previous phases 
are dragged up into regions of larger (mR). This is what we expect given that the 
critical mass in the N —>■ oo limit occurs at mR = 3.203 for all the gapped phases. 
This phenomenon is particularly clear by examination of (mi?) c ^ of the de-confined 
phase for N = 3, 4, 5, 6 in Figure || which shows a steady increase in (mR) c _,d from 
2.6 < (mR) c ^ d < 2.8 for N = 3 to 3.0 < (mR) c ^ d < 3.2 for N = 6. 

It serves as a useful check to compare quantitatively with the results on M 3 x S 1 
in [8]. This can be done by plotting the phase diagram in the (L/R,mL) plane 
and considering the limit L/R — > 0. The phase diagrams for iV = 3,4,5,6 in the 
(L/R,mL) plane are shown in Figure [|. Taking the limit L/R — > in these phase 
diagrams shows precise agreement with the IR 3 x S 1 results in [8] for the values of 
mL at which the transitions occur. Away from the limit of small L/R the phases 
do not take the precise Polyakov loop eigenvalues which are determined for R 3 x S 1 
in [8], rather the eigenvalues spread out slowly around the circle as L/R is increased, 
causing the magnitude |TrP| to decrease. 



4.1 Comparison with lattice results of Cossu and D'Elia 



To connect with strong coupling results it is useful to qualitatively compare the phase 
diagram on S 3 x S 1 to that from the recent lattice simulations in [10]. In [10] the 
authors also consider QCD(Adj) with periodic boundary conditions on fermions for 
the purpose of checking the volume dependence of the phase diagram. They obtain 
results for N = 3 and Nf = 4 (or = 2 Dirac flavours). To compare with their 
results we calculated the phase diagram on S 3 x S 1 in the (mR, R/L) plane as shown 
in Figure [H] (Left). The boundaries of the 2-gap phase correspond to the expected 
lines of transition in the limit of large R, i.e. the M 3 x S 1 result. With permission 
of the kind authors of [10] we show their phase diagram in the (ma, (3) plane in 



Figure |10] (Right), where m is the fermion mass, a is the lattice spacing, and (3 is the 
lattice parameter which goes like the inverse coupling (3 = 2N/g 2 . It is important to 
note that this quantity is frequently confused with L (which is often also called (3), 
however the lattice (3 = 2N/g 2 goes more like 1/L. Obtaining a more quantitative 
relationship between f3 and L would allow for an even better comparison than what 
we have shown. 

In the lattice phase diagram of Figure [TO] (Right) it is unclear whether or not the 
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Figure 10: QCD(Adj) for N = 3, N f = 4: (Left) (mR, L/R). L = 2ttR s i. Only the 
confined phase persists for mR < 3.6; (Right) Results from lattice simulations of Cossu 
and D'Elia [10] on a 16 3 x L c lattice. Here (3 is related to the inverse coupling /? = 2N/g 2 . 



de-confined phase persists into the chiral limit. In the perturbative phase diagram of 
Figure [10] (Left) it does not, however, it is possible that in the case of more strongly 
interacting adjoint fermions the de-confined phase drops down to lower values of mR. 
One thing which may help answer this question is to plot the lattice phase diagram 
in the (mL s , L s j L t ) plane, as it may show a clearer trend. 

It is interesting to compare as well with the lattice phase diagram which results 
from adding a single double trace term to the pure gauge theory action [6,7]. The 
phase diagrams for these theories seem to suggest that the confined phase passes 
through (i.e., the strong and weak-coupling confined phases are connected) when the 
double trace term is strong enough. The double trace term, which is approximately 
the contribution of static massive adjoint fermions, is given by hTr^P for N = 3, 
where larger magnitudes of h correspond to smaller fermion mass. 7 

Another remarkable similarity between the S 3 x S 1 result and the lattice results 
[10] of Cossu and D'Elia is that there is good qualitative agreement for TrP as a 



function of L for a fixed value of the fermion mass. In Figure [II] we show TrP 
as a function of the length of the temporal direction. The S 3 x S 1 result in the 
(L/R, |Proj Zs TrP) plane is in Figure [Ll] (Left). The lattice result of [10] in the 
(Re [TrP] , Im [TrP]) plane for 16 3 x L c lattices with increasing length L c of the 
temporal dimension is given in Figure [11] (Right). What is remarkable about this 
comparison is that the perturbative result appears to capture even some fine details 
of the phase diagram obtained on the lattice including microscopic changes in the 
magnitude of TrP. The sharp discontinuity in TrP at the transition to and from 
the 2-gap phase, and the slow drop in TrP in the de-confined phase as the confined 
phase is approached, seem to agree rather well. 



7 More information on double trace deformations can be found in [31-33]. 



Figure 11: QCD(Adj) for N = 3, Nf = 4: (Left) Results from perturbation theory in 
(L/R, ^Proj^Tri-*) plane. (mR) = 6; (Right) Results from lattice simulations of Cossu 
and D'Elia [10] in the (Re [TrP] , Im [TrP]) plane on 16 3 x L c lattices for (3 = 5.75 and 
am = 0.10. 

It might, at first, seem surprising that agreement between the lattice calculation 
on (S 1 ) and our calculation on S* 3 x S* 1 is so good since these spaces have different first 
homotopy groups. However, 3 of the circles of the torus in the lattice calculation 
are large and so one might expect that the one remaining small circle plays the 
dominant role in determining the phase structure 8 . In addition, it seems that the 
phase diagram doesn't change much as the coupling strength is increased from the 
weakly interacting limit to the strongly interacting limit when considering the phase 
diagram in terms of patterns of Z N breaking. This is not necessarily the case for all 
observables. 



5. Discussion 

It is important to clarify the implications of our results for the issue of volume 
independence. The latter relies on the fact that the theory is in the confining phase 
with unbroken centre symmetry. Our result show that at weak coupling on S* 3 x S* 1 
the confining phase persists in the limit L/R — > as long as the fermion mass is 
below a critical value in units of 1/R, i.e. 1.225/ R for Nf = 1, 3.203/i? for Nf = 2. 
The critical mass increases with increasing Nf. Of course, our result is valid on a 
small S 3 and the interesting question is what happens to the critical mass as one 
moves to strong coupling by taking R > 1/Aq CD . 

The general analysis presented in this paper could just as well be applied to 
many other theories. For example, it is straightforward to consider QCD(Adj) with 
anti-periodic boundary conditions which is of interest since lattice simulations have 

8 We wish to thank Mithat Unsal for useful discussions on this topic. 



been done [34,35], as well as analytical calculations [36], which suggest the presence 
of a conformal window. It would also be interesting to consider symmetric and an- 
tisymmetric representation fermions and compare with related lattice results [37]. 
Interesting also, would be applications to softly broken Af = 4 theory, since there 
are several fermions and scalars and more complicated mass hierarchies are possible 
leading to more complicated phase diagrams. Whether a connection could be made 
with the string theory dual remains to be seen. In particular, it would be interesting 
to understand how the confinement /de- confinement transitions can occur with peri- 
odic boundary conditions for the fermions suggesting that it is not a Hawking-Page 
transition in AdS. The other related issue is how the nature of the transition changes 
as interactions are turned on and one moves to strong coupling. As we mentioned 
in the introduction, for pure gauge theory the 3-loop calculation in [2] shows that 
the transition is first order and occurs at a lower termperature, i.e. larger L/R, than 
the non-interacting Hagedorn transition. Unfortunately it will be very difficult to 
generalize the calculation of [2] to include massive fermions. 

Calculations on S* 3 x S* 1 can be used to define better the extent of perturba- 
tive validity of orientifold planar equivalence [38,39], which is a large N equivalence 
of QCD(Adj) and QCD with symmetric/antisymmetric representation fermions. The 
same can be done for orbifold planar equivalence [40,41], the equivalence of QCD(Adj) 
and QCD with bi-fundamental representation fermions. In both cases a comparison 
of large iV phase diagrams is possible, and perhaps some observables could be com- 
pared. The one- loop effective potential and phase diagrams for QCD(Adj/AS/S) 
on S 3 x S* 1 with massless fermions were computed in the very clearly written pa- 
pers [27,42]. One could also study QCD (fundamental representation fermions) and 
incorporate a finite chemical potential. However, in this case one has to confront the 
sign problem. It would be interesting to compare a phase-quenched result (demand- 
ing a real fermion determinant) on S* 3 x S* 1 with results from the several different 
techniques for dealing with the sign problem in QCD at various coupling strengths, 
a diverse sampling of which can be found in [43-49]. 

There are several things that might be done to allow for better comparison of 
weak-coupling results on S 3 x S* 1 to lattice results. Using the two- loop renormalization 
group equation to give a fitting function for the relationship between the lattice 
parameter j3 = 2N/g 2 and the length L of S* 1 would allow for more quantitative 
comparisons. More lattice results on different volumes and for various iV would 
show if phase diagrams are consistently similar. In particular, we see a pattern 
emerging that suggests that for QCD(Adj) with finite even N there are N/2 gapped 
phases (including the de-confined phase) with the property that |TrP fc | > |TrP'| in 
a k-gap phase for / not a multiple of k, and for N odd there are (N + l)/2 gapped 
phases with this property (narrow regions of additional phases are also possible). 
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On the analytical side interactions might be included in the weak-coupling effective 
potential by working out higher loop orders. In addition, to compare even better with 
lattice results the theory can be put on the torus. The one-loop effective potential 
for QCD(Adj) with massive fermions was computed on R d x T n in [50] (see also [51]). 
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Appendix A: Spherical Harmonics 

In this appendix we collect together results for the spectra of various Laplace 
operators on a sphere S d . 

First of all the scalar Laplace equation is solved by generalized spherical har- 
monics, 



A%(fl) = -e\ s)2 Y £A (h) . 



(A.l) 



The eigenvalues are 



e\ s)2 =£(£ + d- 1)R- 2 , 



(A.2) 




(A.3) 
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Following [52] the eigenvalues of the vector Laplacian on vector fields are 
obtained by separation of the longitudinal (L) and transverse (T) components. The 
spatial gauge field is thus decomposed = Bi + Cj, where Bi is the transverse 
component with Vj-Bj = 0, and Cj is the longitudinal component with C l = V l x- 

The vector Laplacian acting on the longitudinal component Cj has the same 
spectrum as for the scalar Laplacian because 

(V*V4 - R j k ) V jX = V fe (V l V lX ) , (A.4) 

except that the £ = mode is missing: 

e^ L)2 = £(£ + d- l)iT 2 , (A.5) 

for £=1,2,..., and 



{V , L) _ (2£ + d-l)(£ + d-2)\ 



The eigenvalues and degeneracy of the vector Laplacian on the transverse com- 
ponents Bi are 

e^ 2 = {£{£ + d - 1) + d - 2) IT 2 , (A.7) 

and 

,(,,T) l(l + d-\){2t + d-W + d-Z)\ 

d * = (d-2)\(£+l)\ • (A ' 8) 

for £= 1,2,.... 

The Laplacian on spinors is given by 

A (/) = TV'ViVj- = g ij V i V j - \R , (A.9) 
where 1Z is the scalar curvature of the sphere and 

V i = 9 i + r i (A. 10) 

is the covariant derivative on spinors with spin connection IV The eigenvalues and 
degeneracy [53-55] are 

e ^ 2 = (€ + |) 2 ir 2 , (A.ii) 

and 

,(/)_ 2(rf + £-2)! 



d ' ~ (e-i)\(d-i)\ ' (A ' 12) 



where £=1,2, 
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Appendix B: The Abel-Plana Formula 



In this appendix, we prove the version of the Abel-Plana formula quoted in the 
main text ( 2.16|) . The idea is to represent the sum on left-hand side as a contour 
integral: 



e=o 



(B.l) 



where C is the contour illustrated in Fig. [L2[ In our case, the function f(z) has 
square root branch points at z — ±imR. Using the above we have 



00 PC 

£^+§)- / 



dx f(x) 



2i 



oo+ie 



dz f(z)(tan(nz) — i) — 



c.c. 



(B.2) 



Now we can rotate the contour here that runs from the origin out to infinity over the 
poles so that it runs from the origin to ioo to the right of the branch point of f(z) 
at imR. This gives the right-hand side as 



' f(ix + e) 
— 1 1 dx — ; C.C. 



+ 1 



(B.3) 



Hence, we have proved that 



£=0 



+ 1 



(B.4) 



which is the formula (|2.16|) in the text. 



v 



Figure 12: The contour used in the derivation of the Abel-Plana formula. 
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Appendix C: Numerical minimization 

In the finite N calculations of the phase diagram it was necessary to numeri- 
cally minimize the effective potential. This was not trivial as it is difficult to find 
an algorithm that will always find the global minimum of an arbitrary function. 
Two techniques which proved most useful were the Random Search and Differential 
Evolution numerical minimization routines implemented in Mathematica, which are 
reviewed in [56]. Increasing the number of search points improves the chances of 
obtaining the global minimum. Somewhat surprisingly in most cases just 5 search 
points were enough to obtain the minimum accurate to around 13 digits for N — 3. 
It is only slightly less accurate when considering N = 4, 5, 6. The addition of the 
mass term to the effective potential actually makes the minimization easier. The 
pure gauge theory plots required many more search points, between 500 and 1500. 

It was not possible to obtain reliable minimization of the effective potential when 
there were terms represented as an infinite series. Therefore it was necessary to put 
all infinite series in a non-series expression, either by solving them, or converting 
them into integral forms using the Abel-Plana formula. 

As this type of calculation doesn't have much history (however, the same proce- 
dure was used in [8]) there is a lot of room for improvement in technique. 

It is additionally important to perform checks of the saddle point approximation 
by plotting the relevant observables as a function of the configuration space of the 
9i as discussed in the finite N section. 
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